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In certain situations the state of a quantum system, after transmission through a quantum channel, 
can be perfectly restored. This can be done by "coding" the state space of the system before 
transmission into a "protected" part of a larger state space, and by applying a proper "decoding" 
map afterwards. By a version of the Heisenberg Principle, which we prove, such a protected space 
must be "dark" in the sense that no information leaks out during the transmission. We explain 
the role of the Knill-Laflamme condition in relation to protection and darkness, and we analyze 
several degrees of protection, whether related to error correction, or to state restauration after a 
measurement. Recent results on higher rank numerical ranges of operators are used to construct 
examples. In particular, dark spaces are constructed for any map of rank 2, for a biased permutations 
channel and for certain separable maps acting on multipartite systems. Furthermore, error correction 
subspaces are provided for a class of tri-unitary noise models. 



I. INTRODUCTION 

We consider a quantum channel of finite dimension through which a quantum system in some state is sent. The 
output consists of another quantum state, and possibly some classical information. We are interested in the question 
to what extent the original quantum state can be recovered from that state and that information. In particular, we 
investigate if there are subspaces of the Hilbert space of the original system, on which the state can be perfectly 
restored. 

In the literature a hierarchy of such spaces, which we shall call protected subspaces here, has been described. The 
strongest protection possible is provided in the case of a "decoherence free subspace" [1-4] . In this case the channel 
acts on the subspace as a isometric transformation. All we have to do in order to recover the state, is to rotate it 
back. 

The next strongest form of protection occurs when the channel acts on the subspace as a random choice between 
isometrics, whose image spaces are mutually orthogonal. Then by measuring along a suitable partition of the output 
Hilbert space, it can be inferred from the output state which isometry has occurred, so that it can be rotated back. 
This situation is characterized by the well-known Knill-Laflamme criterion, [5, 6] and the protected subspace in this 
case is usually called an error correction subspace. 

The weakest form of protection is provided in yet a third situation, which was encountered in the context of quantum 
trajectories and the purification tendency of states along these paths [7] . In this case the deformation of the state is 
not caused by some given external device, but by the experimenter himself, who is performing a Kraus measurement 
[8]. Also in this case the "channel" acts as a random isometry, but the image spaces need not be orthogonal. It is now 
the measurement outcome (not the output state), that betrays to the experimenter which isometry has taken place. 
Using this information, he is able to undo the deformation of the component of the state that lies in the subspace 
considered. 

It should be emphasized that the latter form of protection is far from a general error correction procedure. The 
experimenter only repairs the damage that he himself has incurred by his measurement. 

Nevertheless, the above situations seem mathematically sufficiently similar to deserve study under a common title. 

In all these three cases the experimenter learns nothing during the recovery operation about the component of the 
state inside our subspace. In this sense these subspaces can be considered "dark" , and this darkness is essential for 
the protection of information. Our main result (Theorem 3) is concerned with the equivalence between protection 
and darkness, which is a consequence of Heisenberg's principle that no information on an unknown quantum state 
can be obtained without disturbing it (Corollary 2). 

The question arises, for what channels protected subspaces are to be be found. We consider several examples in 
their Kraus decompositions. In each decomposition, we look for subspaces on which the channel acts as a multiple 
of an isometry, to be called a homometry here. Obviously, every (Kraus) operator A acts homometrically on a one- 
dimensional space (D'0; its image CAtp is another one-dimensional space, and the shrinking factor is ( Aip, Aip ) = 
However, one-dimensional spaces are useless as coding spaces for quantum states. What we shall need. 
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therefore, is the recent theory of higher rank numerical ranges [9, 10]. With the help of this we shall be able to 

construct several examples. 

The paper is organized as follows. A brief review of basic concepts including channels and instruments is presented 
in section II. We discuss Heisenberg's principle in Section III. and prove our main Theorem, Theorem 3 in Section IV. 
In subsequent sections we analyze different forms of protected subspaces and compare their properties. In section V 
we review the notion of higher rank numerical range and quote some results on existence in the algebraic compression 
problem. Some examples of dark subspaces are presented in section VI, while an exemplary problem of finding an 
error correction code for a specific model of tri-unitary noise acting on a 3 x ii' system is solved in section VII. 



II. CHANNELS AND INSTRUMENTS 



Let H be a finite-dimensional complex Hilbert space, and let B{H) denote the space of all linear operators on Ti. 
We consider H as the space of pure states of some quantum system. By a quantum operation or channel on this 
system we mean a completely positive map $ : BCH) B{'H) mapping the identity operator 1 = I-h to itself. The 
map $ describes the operation "in the Heisenberg picture", i.e. as an action on observables. Its description "in the 
Schrodinger picture", i.e. as an action on density matrices p, is described by its adjoint The maps $ and $* are 
related by 

V,VxeB(w) : tr($*(p)X) = tr(p$(X)) . 
We note that the property $(1) = 1, which we require for is equivalent to trace preservation by $*: 

tr($*(p)) = tr($*(p) • 1) = tr(p • $(1)) = tr(p • 1) = tr(p) . 
By Stinespring's theorem, every channel $ : B{Ti) BiTi) can be written as 

^{X) = V\X(S)1m)V , (2.1) 

where V is an isometry H ^ H. ® M. iov some auxiliary Hilbert space M.. The minimal dimension r oi M. admitting 
such a representation is called the Choi rank [11, 12] of 

Any Stinespring representation of $ naturally leads to a wider quantum operation 

* : B{H) B{M) B{H) -.X^Y^V^X® Y)V , (2.2) 

which can be interpreted (in the Heisenberg picture) as the result of coupling the system to some ancilla having 
Hilbert space A4. 

Thus Stinespring's representation (2.1) can be symbolically rendered as in Fig. 1. 




B{M) 

FIG. 1: Stinespring's dilation of $ seen as coupling to an ancilla M 



In this picture, the cross stands for the substitution of 1m (in tiie Heisenberg picture, reading from right to left), or 
the partial trace (in the Schrodinger picture, reading from left to right). Physically, it corresponds to throwing away, 
or just ignoring, the ancilla after the interaction. In the picture, the fact that is a compression, i.e. = V"'' • F for 
some isometry V, is symbolized by the triangular form of its box. 
Now, by blocking the other exit in Fig. 1, we obtain the conjugate channel [13], 

^'^ : B{M) B{n) : Y ^ *(!« O F) = V^{ln O Y)V . 

See also Fig. 2. 

The main message of this paper is the following. The conjugate channel can be viewed as the flow of information into 
the environment. By Heisenberg's Principle, to be explained below, such a flow prohibits the faithful transmission of 
information through the original channel In particular, if the information encoded in some subspace of H is to be 
transmitted faithfully, nothing of it is visible from the outside: protection implies darkness. The degree of protection 
(decoherence free, strong or weak) is related to the degree of darkness, for which we shall define some terminology. 



FIG. 2: The conjugate channel . 



B{M) 



Any orthonormal basis / = (/i, . . . , /„) in A4 corresponds to a possible von Neumann measurement IVj.- on the 
ancilla, which maps a density matrix p on to a probabihty distribution {{ fi, pfi) ,{ f2, pf2) ,-■•,{ frm pfm)) 
on {1,2, ... ,m}. (Cf. Fig. 3.) In the Heisenberg picture this is the map from the algebra (D™ with generators 
ei = (1,0,0,...,0), 62 = (0,1,0,...,0), e„ = (0,0,...,0,1), to B{M), given by 

11/ : \fi){fi\ ■ 




FIG. 3: Von Neumann measurement on M. 

In FIG. 3 the abelian algebra (D™ is indicated by a straight line since it only carries classical information. Quantum 
information is designated by a wavy line. 
Let us now denote by // the "partial inner product map" 

n^M^H: ip(E)e^ {f,0)ip , 

and let us write 

Ai:=If,V eB{n). 

Then since lj,XIf- = X (g) \fi){fj\, we obtain a decomposition of $ along the basis as follows: 

m mm 

$(x) = nx<^iM) = ® = = E^i^^j- • (2-3) 

i—1 i—1 i—1 

This is a Kraus decomposition of $. Combining the coupling to the ancilla with a von Neumann measurement on the 
latter, we obtain an instrument in the language of Davies and Lewis [14]: 

*/ :fi(H)«)C'"^fi(7W) : X ei V\X ^ \fi}{fi\)V = AjXAi . (2.4) 
The isometric property of V is now expressed as 

m 

V^V = J2aIA, = 1 . (2.5) 

III. HEISENBERG'S PRINCIPLE OR OBSERVER EFFECT 

In quantum mechanics observables are represented as self-adjoint operators on a Hilbert space. When A and B 

arc commuting operators, then they possess a common complete orthonormal set of eigenvectors. Each of these 
eigenvectors -0 determines a state which associates sharply determined values to both observables A and B. 

But when A and B do not commute, such states may not exist. This important property of quantum mechanics was 
first discussed by Heisenberg [15], and is called the Heisenberg Uncertainty Principle. It was formulated by Robertson 
[16] in the form 



a^{A) ■ a^{B) > i] ( V, {AB - BA)^;) \ 
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Here cr^(X) is the standard deviation of X in the distribution induced by tfj. Already in the very same paper, 
Heisenberg introduced a second and very different principle, which is sometimes designated as the "Observer Effect" , 
and which we shall call the Heisenberg Principle here. Roughly speaking, it says that: 

if A and B do not commute, 

a measurement of B perturbs the probability distribution of A. (3.1) 

In the first half century of quantum mechanics, physicists, including Heisenberg himself, were satisfied with this 
formulation, and even considered it more or less identical to the Uncertainty Principle above. 

In recent years it was realized that in fact we have here two different principles. Good quantitative formulations 
have been given of the Heisenberg Principle (for example [17, 18]). For the purpose of the present paper we are 
satisfied with a qualitative ('yes-or-no') version. 

Let us first note that the formulation of the principle needs sharpening. As it stands, the condition is not needed: 
already in the trivial case that A ^ B measurement of B changes the probability distribution of A. Indeed changing 
the probability distribution of an observable is the very purpose of measurement! And also, when A and B commute, 
but are correlated, then gaining information on B typically changes the distribution of A. A characteristic property of 
quantum theory only arises if we require that the outcome of the measurement of A is not used in the determination 
of the new probability distribution of B. Even then, some states may go through unchanged. 

Corrected for these observations, the Heisenberg Principle reads: 

For noncommuting A and B we cannot avoid that, 

for some initial states, a measurement of B changes the distribution of A, 

even if we ignore the outcome of the measurement. (3-2) 

The contraposition of the statement turns out to be mathematically more tractible: 

If the probability distribution of A is not altered in any initial state 
— by us performing some measurement and ignoring its outcome — 

then the object measured must commute with A. (3.3) 

In this form it is sometimes called the 'nondemolition principle'. 

Now let us make this statement precise. We start with a self-adjoint operator A on Ti. Its distribution in the state 
p is determined by the numbers tT:{pg{A)) when g runs through the functions on the spectrum of A. Then some 
quantum operation is performed which on BiJ-L) is described by a completely positive unit preserving map We 
require that for all states p and all functions / 

tv{^*{p)g{A))=tv{pg{A)) , 

which is equivalent to 

^{g{A)) = g{A) . 

I.e.: all elements of the *-algebra A consisting of functions of A are left invariant by Let us denote the commutant 
of A by A', 

A' = {XG B{n) I WreA ■ XY = YX} . (3.4) 
Now, the quantum operation $ is due to a measurement, so it is actually of the form 

$(x) = e(X(g)i), 

where © : B(W)(g)(D'" — > B(W) is some instrument whose outcomes, labeled 1, 2, . . . , m, in the state p have probabilities 
Pi,P2,--- ,Pm to occur, where 

Pj = tr(pe(l Cj)) , 

and where tr[pQ{X ej))/Pj is the expectation of X, conditioned on the outcome j. (This situation is comparable 
to, but more general than, that of 5"^ in (2.4).) Here (D™ is the algebra of measurement outcomes. Generalizing to 
arbitrary A, we may now formulate the Heisenberg Principle as follows. 
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Proposition 1 (Heisenberg Principle.) LetTL he a finite dimensional Hilbert space, andB some finite dimensional 
*-algebra. Let A be a sub-*-algebra of B{'H), and let Q be a completely positive unit preserving map B{T-L)®B B{l-L). 
Suppose that for all A & A we have 

e(A0 i) = A. 

Then for all B & B 

e(l ®B)&A! . 

Proof: For any density matrix p on H, define the quadratic form Dp on B by 

Dp{X,Y) := tvp{Q{X*Y) - Q{Xy@{Y)) . 

By the Cauchy-Schwartz inequality for the completely positive map G this quadratic form is positive semidefinite. 
By assumption we have for all A & A: 

Dp{A(^l,A®l) = trp(e(A*A(g) 1) - 9(^0 l)*e(A(g) 1)) 
= trp{A*A O 1 - (A O O 1)) = . 

It then follows from the Cauchy-Schwartz inequality for Dp itself that Dp{A 1, 1 (g) -B) = for all B G B. But then 

trp{Ae{l (g) B)) = trp{e{A 1)6(1 ® B)) = trp(e((.4 ® 1)(1 ® B))) = trp(e((l (g) B){A (g) 1))) 
= trp(e(l (g B)e{A (g 1)) = tr/3(e(l g) B)A) . 

Since this holds for all p, it follows that Q{1 g) B) commutes with A. □ 

By taking A and B abelian, say A generated by some observable A, and B = (D™ as above, and by choosing for 6 
some instrument giving information about B, wc obtain a statement of the type (3.3). 

But there are other possible conclusions. We may choose A = B{H), so that A' = <C ■ 1-h- Then the Heisenberg 
principle says that, if we wish to make sure that any possible state p on W be unchanged by our measurement, no 
information at all concerning p can be gained. This is expressed by the following corollary and FIG. 4. 

Corollary 2 In the situation of Proposition 1, if for all A e BiTi) we have 

e(A(g) 1) = A , 

then there is a positive normalized linear form a on B such that for all B £ B: 

9(1 g)S) =a{B) ■ In ■ 
Indeed, the expectation of an outcome observable, 

tr(e»(l B) = tr(pe(l B)) = ir:{plu) ■ tv{aB) = tr(aB) 
does not depend on p (see FIG. 4.) 



IV. PROTECTION AND DARKNESS: THE KNILL-LAFLAMME CONDITION 

Let £ be a complex Hilbert space of dimension smaller than that of 7i, and let C : £ — > W be some isometry. The 
range of C is a subspace of H, isomorphic with C. Let T : B{H) — » B{£) denote the compression map 

r{x) = c^xc . 

Note that T is completely positive and identity-preserving. Compression maps are a convenient way of describing 
subspaces of a Hilbert space in the language of operations. Note that the operation F* (in the Schodinger picture) 
embeds density matrices on C into the range of C: 



r*{p) = CpC^ 



FIG. 4: Heisenberg's Principle as an implication between diagrams 




Physically, T is to be viewed as the "coding" operation. 

Definition. Wc say that F (or the subspacc CC of H) is protected against a channel $ : B{'H) — > B{'H) if F o $ is 
right-invertible, i.e. if there exists a "decoding" operation A : B{£) BiTL) such that 

Fo3>oA = ide(£) . (4.1) 

By virtue of (2.1) we may picture this state of affairs as in Fig. 5. 

The subspace will be called weakly protected against an instrument ^ / : BiTi) (g) (D™ BiTL) liVo'^ f'ls right-invertible, 
i.e. if there exists a decoding operation A/ : B{C) — > B(W) (g) C" such that 

Fo^-joA/ =idB(£) . (4.2) 

This is symbolically rendered in Fig. 6. The difference with Fig. 5 is that, in the case of weak protection, it is 
allowed to use the measurement outcome in the decoding. In the figure the classical information consisting of the 
measurement outcome, is symbolized by a straight line. 




FIG. 6: Weak protection of F against ^ / 



The above notions concern protection of information. Now we consider its availability to the external world. 

Definition. Let 5*/ : B{TL) ® (D™ B{H) denote a quantum measurement (instrument) as described in (2.4). 
The subspace C£ c H (or the compression operation T = ■ C), will be called dark with respect to if for all 
i = 1, . . . ,m we have 

ro*/(l(g)ei) e C- l£ . (4.3) 
This condition can be written in an equivalent form, 

CUJa,C = Xirlc for z = 1, . . . , m . (4.4) 

The subspace C£ will be called completely dark for a channel $ : B{H) BiTL) if it is dark for all Kraus measurements 
^ f obtained by choosing different orthonormal bases in the ancilla space of some Stinespring dilation of i.e. 

"^YeBiM)- ro^'(loF) e C- l£ . (4.5) 

In terms of Kraus operators this is equivalent with the Kmll-La.fia.mme condition: 

C^'AIAjC = aij -Ic for j = 1, . . . , m . (4.6) 
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Interpretation: From (4.3) and (4.4) we see that, if the von Neumann measurement along / is performed, the 
measurement outcome i has the same probabihty p(r o ^1/^(1 (g) Cj)) = p{C''^ A^AiC) = Xi, in all system states p, i.e. 
no information concerning the state p can be read off from the /-measurement on the ancilla. 

Complete darkness (i.e. (4.5) or the equivalent Knill-Laflanime condition (4.6)) says that no information whatsoever 
concerning the input state reaches the ancilla. Mathematically, the Knill-Laflamme condition says that the range of 
the conjugate channel lies entirely in the center (D • l£ of B{C). Let us emphasize again that if the space C satisfies 
the conditions (4.6) for a map ^ represented by a particular set of the Kraus operators {Ai}^i, then C also satisfies 
them for any other set of Kraus operators {-B^}™ i, used to represent the same map VP. 

Note also that the set of conditions (4.6), which express complete darkness, naturally defines a state a, on the ancilla 
by a relation 

To^{l(g)Y) =tT{aY)-lc . (4.7) 

satisfied by any Y . This quantum state acting on an auxiliary system is called the error correction matrix, since the 
density matrix aij appears in eq. (4.6). Observe that the density operator a depends only on the map and not 
on the concrete form of the Kraus operators Ai, which represent the map and determine the matrix representation 
aij of a. Relations between matrix elements of the same state represented in two different basis are governed by the 
Schrodinger lemma [12], also called GHJW lemma [19, 20]. 

We are now going to prove the equivalence of protection and darkness. In the case of strong protection and complete 
darkness this reproduces and puts into perspective the result of Knill and Laflamme [6] In that case, if the state a 
is pure, then the decoding operation A can be realized by a unitary evolution, Hence the purity constraint for the 
error correction matrix, a = a^, is the correct condition for a decoherence free subspace [21] - see also the proof 
of Theorem 3. As a quantitative measure, which characterizes to what extent a given protected space is close to a 
decoherence free space, one can use the von Neumann entropy of this state, S = — Tra In a. This code entropy [22] 
is equal to zero if the protected space is decoherence free or if the information lost can be recovered by a reversible 
unitary operation. Observe that the code entropy S characterizes the map ^E* and the code space C, but does not 
depend on the particular Kraus form used to represent 5". 

In this way we have determined a hierarchy in the set of protected spaces. Every decoherence free subspace belongs 
to the class of completely dark subspaces, which correspond to error correction codes. In turn the completely dark 
subspaces form a subset of the set of dark subspaces - see Fig. 7. 



decoherence 
free ^ 



protected subspaces 

FIG. 7: Sketch of the hierarchy of protected subspaces. 




Theorem 3 (Equivalence of Protection and Darkness) Let Ti., A4, and C be finite dimensional Hilhert spaces. 
Let C : C^H and V -.H^H^M be isometrics, and let <&, ^E" and ^/ be as defined in (2.1), (2.2) and (2.4). Then 
CC is weakly protected against the instrument ^ f if and only if CC is dark for . It is strongly protected against $ 
if and only if it is completely dark for <&. 

Proof: 

First assume that CC is strongly protected against <f>, i.e. (4.1) holds for some decoding operation A. Let ^{X) — 
^'(X (g) 1) for some compression ^. Define 

e : BiC) (E) BiM) B{C) -.X^Y^To *(A(X) ® Y) . 

Then Q{X ®1)^X for aU X e B{C), and by Corollary 2, since A(l) = 1, 



r o *(i (g) y) = 9(1 ® r) e c • 1 
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so (4.5) holds, and CC is completely dark for 

Conversely, suppose that C£ is completely dark for and let a denote the density matrix given by (4.7) Then we 
may diagonalize: 

m 

tv{aY)=^ai{fi,Yfi) 

i=l 

for some orthonormal set (/i)^i (with m' < m) of B{Ai) and positive numbers oi, 02, . . . , a^,' summing up to 1. Now 
let Ai := If.V. Then for all ip & C: 

{AiCiP,AjCiP) = {lf,VCi^,If^VCij) 

= (^,Ctl/t(i^|/.)(/.|)yc'V) 

= aiSij ■ \\tpf . 

So the ranges of AiC and AjC are orthogonal for i ^ j and Ai is homometric on CC Now define Dj for i = 1, 2, . . . , m' 
on these orthogonal ranges by 

Di(p = if _L Range (j4jC), DiAiCip = y/oitp . 
{Di "rotates back" the action of AiC.) Let A denote the operation 

m' I 7n' \ 

A(Z) := ^ D\ZDi + p{Z) \1h-Y. ■ 
i=i \ j=i J 

for some arbitrary state p on B{£). (The term with p is intended to ensure that A{lc) = Then we have for all 
Z e B(£): 

m' m' 

ro$oA(Z) = Y^Y^C'^AXolZDiAjC 

j=i i=i 

= —C''A]AiCZC^AlAjC = J2 ^ijo-i^ = ^ ■ 

j=l i=l * ij=l 

So CC is strongly protected against <1> by (4.1). 

Now let us prove the equivalence between weak protection and darkness. Assume that C£ is weakly protected 
against i.e. (4.2) holds for some A/ : B{C) B{n) O C™, say A/(X) = Ejli^/l^) ® ^i- Define 6 : 
B{n) (D™ ^ B{n) by 

m 

e(X ® g) := ^5(j)r o 4'/(A^ (X) ® e,) . 

Then by (4.2), Q{X ®1) = X for all X € B{C). Hence by Corollary 2, 

To-^f{l®ei) = Q{l®ei) G B{n)' = C • 1^ . 

So (4.3) holds, and CC is dark for 

Conversely, assuming that CC is dark for ^f/, then AiC is homometric on C by (4.4), and we may define Di : H ^ C 

by 

DiAiCi) := ^/Xii) iftpeC, Dup = if ± Range (AC) . 
(Briefly: A = CU^/VAl if A; ^ 0, zero otherwise.) Define the decoding operation A/ : B{C) B{n) O by 

m 

^/(^) {DjZDi + {In - D\Di)p{Z)) 
1=1 
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for some arbitrary state p on B{C). Then, for Z e B{C): 

/ m N 

ro*;oA/(z) = ro*W ^(d|za + (i--d/aM2)) oej 



1=1 



= C^y^ J^iDjZDi + (1 - DjDi)p{Z)) ® \fi){fi\ VC 



m m m 

^Ct^DjZAAiC = ^ -{C^AjAiC)Z{C^AjAiC) = Y,hZ = Z . 



□ 



V. COMPRESSION PROBLEMS AND GENERALIZED NUMERICAL RANGE 

For a given channel $ : — > we are interested in the protected subspaces of Ji. These are the subspaces 

on which the compressions of aIAj act as scalars. In this section wc review this compression problem. 

Let T be an operator acting on a Hilbert space Ti. of dimension n, say. For any fc > 1, define the rank-k numerical 

range of T to be the subset of the complex plane given by 

Afe(T) = {A e (D : C^TC = Al for some C : (D*^ ^ U), (5.1) 

The elements of A;; (T) can be called "compression- values" for T, as they are obtained through compressions of T to 
a fc-dimensional compression subspace. The case k = 1 yields the standard numerical range for operators [23] 

Ai(T) = {(V|TV) : IV) e W , im = 1}. (5.2) 

It is clear that 

Ai(T)DA2(T)D...D A„(T). (5.3) 

The sets Afc(T), fc > 1, are called higher-rank numerical ranges [9, 24]. For any normal operator acting on 7i„ 
this is a compact subset of the complex plane. For unitary operators this set is included inside every convex hull 
(coF), where F is an arbitrary (n + 1 — fc)-point subset (counting multiplicities) of the spectrum of T [9]. It was 
recently shown that for any normal operator the sets Afc(T) are convex [25, 26] while further properties of higher rank 

numerical range were investigated in [27-29]. 

The higher rank numerical range is easy to find for any Hermitian operator, T = acting on an n-dimensional 
Hilbert space H. Let us quote here a useful result proved in [9]. 

Lemma 4 Let xi < X2 < ■ ■ ■ < Xn denote the ordered spectrum (counting multiplicities) of a hermitian operator T. 
The rank-k numerical range of T is given by the interval 

Afc(r) = [xk,Xn+l-k] , (5.4) 

Note that the higher rank numerical range of a hermitian T is nonempty for any k < int[(T?,+ 1)/2] . Let us demonstrate 
an explicit construction of a compression to (D^ which solves equation (5.1) for a Hermitian matrix T of size n — A. 
The latter's eigenvalue equation reads T\(j)i) = Xi\(pi). Choose any real A € ^2{T) = [x2,X3]. It may be represented 
as a convex combination of two pairs of eigenvalues {xijX^} and {x2, x^} - see Fig. 8a. Writing 

A = (1 — a)xi + axs = (1 — b)x2 + bx^ (5.5) 

one obtains the weights 

•^"^1 -2/) A u -2/) a\ 

a = =: sm C'l and o = =: sm 02 (5.6) 

.T3 — Xl X4 — X2 

which determine real phases 9i and 62- These phases allow us to define an isometry C : (D^ ^ H by 

„ J ei 1-^ cos6li|(/>i) + sin6li|(^3) , , 

{€2 ^ COS02|'^2)+sin02|04) ' ^ 
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Observe that 

(ei,C'lTCei ) = cos6'ia;i(0i|?/'i) + sin6'ia;3(03|?/'i) = (1 - a)xi+ax3 = A. (5.8) 

Similarly, we have {e2,C'<TCe2) = A. Further, we also have ( ei, C^TCea ) ( 62, CtrCei ). It follows that 
C^TC = 1, and the isometry (5.7) provides a solution of the compression problem (5.1) as claimed. Note that one 
can select another pairing of eigenvalues, and the choice {xi,X4^} and {x2,X3} allows us to get in this way another 
subspace C'C spanned by vectors obtained by a superposition of states with |04) and |02) with jc^g) respectively. 



a) b) 




FIG. 8: Standard numerical range Ai and higher rank numerical range A2 for a) Hermitian operator T of size 4 and b) 
non-degenerate unitary U € C/(4). Observe similarity in finding the weights a and b used to construct superposition of states 
forming the subspace CC in both problems. 

For a given operator T one may try to solve its compression equation (5.1) and look for its numerical range Ak{T). 
Alternatively, one may be interested in the following simple compression problem: For a given operator T find all 
possible subspaces CC of a fixed size k which satisfy (5.1). 

Furthermore, it is natural to raise a more general, joint compression problem of order M. For a given set of 
M operators {Ti, . . .T^} acting on 7i„ find a subspace C£ of dimensionality k which solves simultaneously M 
compression problems: 

C^T,nC = X,n 1 for m = 1, . . . , M . (5.9) 
Note that all compression constants. Am G Ak{Tm), can be different, but the isometry C needs to be the same. 

VI. DARK SUBSPACES 

In this section we provide several results concerning existence of darks spaces for several classes of quantum maps. 

A. Random external fields 

Consider a noisy channel ^ given by 

r 

<^u{X)=Y,qMxU,, (6.1) 

1=1 

where all operators Ui are unitary while positive weights qi sum up to unity. Such maps are called random external 
fields [30] or random unitary channels. The standard Kraus form (2.3) is obtained by setting Ai — ^/qlUi. 

In this Kraus decomposition the whole space, and hence every subspace, is dark. This corresponds to the fact that 
the choice between the unitaries, which is made with the probability distribution (gi, . . . , g^), gives no information on 
the quantum state. And indeed, knowledge of the "external field", i.e. of the outcome «, permits us to undo, by the 
inverse of Ui, the action of the channel. 
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B. Rank two quantum channels 

Let us now analyze a rank two channel, 

p' = $2(p) = A^pAl + A2pA\ , (6.2) 

Lemma 5 For any Kraus representation of any rank-two channel acting on a system of size N there exist a dark 
subspace of dimension k = int[(A'' + l)/2]. 

Proof. Wc need to solve a joint compression problem (5.9) of order two, for two Hermitian operators Ti = A\Ai 

and T2 = A2A2. Due to Lemma 4 there exists a subspace Pk of dimension k = int[(A^ + l)/2] which solves the 
compression problem for the Hermitian operator Ti of size N. It is also a solution of the compression problem for the 
other operator, since the trace preserving condition implies T2 = 1 — Ti . □ 

C. Biased permutation channel 

Consider a quantum map acting on a system of arbitrary size n described by the Kraus form (2.3). Let us assume 
that all Kraus operators are given by 'biased permutations' 

Ai = Pi^/D'i, i = l,...,r. (6.3) 

where Di is a diagonal matrix containing non-negative entries, and Pi denotes an arbitrary permutation of the A''- 
element set. Hence all elements of the POVM form diagonal matrices, 

Ti = A\Ai = ^iPlPi^/Wi = Di , (6.4) 

in general not proportional to identity. Note that the Kraus operators defined in this way need not to be Hermitian. 
To satisfy the trace preserving condition (2.5) we need to assume that X^^^i -Dj = 1. Let us define an auxiliary 
rectangular matrix of size r x N, namely Sim '■= {Di)mm > 0. Then the above constraints for the matrices Di is 
equivalent to the statement that S is stochastic, since the sum of all elements in each column is equal to 1, 

r 

Y,Sim = l for m=l,...,N. (6.5) 

i=l 

A map described by Kraus operators fulfilling relations (6.3) and (6.5) will be called a biased permutation channel. 

We are going to construct a dark space for a wide class of such channels. For simplicity assume that the size of the 
system is even, N = 2k. Let us additionally assume that all elements in each row of B are ordered (increasingly or 
decreasingly) and that the matrix S enjoys a symmetry relation, 

Si,m + Si^n-m+i = const =: Ai for i = l,...,r; m = l,...,k = n/2 . (6.6) 
Then the numbers Aj can be defined by a sum of the entries in each row, Aj = S^m=i ^i'^- 

Lemma 6 Assume that a biased permutation channel acting on a system of size N = 2k possesses the symmetry 
relation (6.6). Then it has a dark space of dimension k = n/2. 

Proof. We need to find a joint compression subspace for the set of r elements of POVM given by diagonal 
matrices Di, with i = 1, . . . ,r. Since these matrices commute, they have the same set of eigenvectors, denoted by 
\vm)-i m = 1, . . . ,N. Due to symmetry relation (6.6) we know that the barycenter of each spectrum, A, belongs to 
the higher rank numerical range, Ak(Di). Furthermore, this relation shows that (for any i) the number Aj can be 
represented as a sum of two eigenvalues oi Di with the same weights, A^ = ^{Di)mm + ^{Di)m'm' with m' = n+l — m. 
By construction this property holds for all operators Di, i = 1, . . .r. Hence the general construction of the higher 
order numerical range for Hermitian operators [10] implies that the subspace 

k ^ 
Ck ■■= where IV'j) := — + |^^i-i+iv)) (6.7) 



12 



fulfills the joint compression problem for all operators Ti = Di, i = 1, . . 

a 



. r. Hence this subspace is dark as advertised. 



To watch the above construction in action consider a three biased permutation channel acting on a two qubit 
system. Hence we set r = 3 and N = 4, and assume that five real weights satisfy < a < b < x/2 < 1/2 and 
< c < d < x/2. They can be used to define the channel by a stochastic matrix S 



b X — b X — 
S = \ c d X — d X — c 
,a' b' b" a" 



(6., 



where a' = 1 — a — c, b' = 1 — b — d, a" = 1 — 2x + a + c and b" = 1 — 2x + b + d. Note that this matrix satisfies the 
symmetry condition (6.6), the elements in each row are ordered, while mean weights in each row read Ai = A2 = x/2 
and A3 = 2(1 — x). 

To complete the definition of the channel we need to specify three permutation matrice of size four. For instance let 



us choose Pi = P(i,2,3,4), P2 = P(i. 



2), (3,4) 



and Pa = P(i 



(1,4,3,2) 



where according to the standard notion, the subscripts 



contain the permutation cycles. Then the biased permutation channel is defined by the three Kraus operators 



f Vb 



y/x^ 

\^/a / 



, A2 



f Vd 






\ 




\/a; — c 
\0 y/x-d J 



, A3 



( 
































V 





^/V' 


/ 



(6.9) 



which satisfy the trace preserving condition (2.5). 

Since the barycenter A; of the spectrum of the POVM element Ti = Di (given by a row of matrix (6.8)), is placed 
symmetrically, in all three cases it can be represented by a convex combination of pairs of eigenvalues with weights 
equal to 1/2. Thus we define two pure states 



|Vi):=^(K) + K)) 



IV2) :=-^(h) + h)), 



(6.10) 



and the two dimensional subspace spanned by them, C = + |V'2)(V'2|- It is easy to verify that the subspace 

C satisfies C^TiC = Ail = C^T2C while C^'T^C = A3I so this space is dark. Note that the subspace CC cannot be 
used to design an error correcting code since C^^j A2C ^ (D • 1. 



D. Composed systems and separable channels 

Consider a bipartite system of size n = ua x ns- A quantum operation $ acting on this bipartite system is called 
JocaJ, if it has a tensor product structure, $ = <I>^ (E" where both maps $^ and $b are completely positive 
and preserve the identity. If for both individual operations, $^ and $b, there exist protected subspaces Ck and Qi 
respectively, then the product subspace Ck ® Qi of size kl is also a protected subspace for the composite map $a ^^^b- 

Similar protected subspaces of the product form can be constructed for a wider class of separable maps (see e.g. 
[12]), 

r 

p' = <S>*{p) = J^i^i ® Bi)p{Ai ® Bi)l (6.11) 

i=l 

Assume that a subspace Ck & Hna is a solution of the joint compression problem for the set of r operators 
while a subspace Qi G Hnb does the job for the set of r operators -BjSj. It is then easy to see that the product 
subspace Ck fX" Qi of dimension kl is a dark subspace for the separable map (6.11). 

It is straightforward to extend lemmas 3 and 4 for separable maps acting on composite systems and apply them 
to construct protected subspaces with a product structure. On the other hand, if for certain problems such product 
code subspace do not exist, one may still find a code subspace spanned by entangled states. Such a problem for the 
tri-unitary model is solved in following section. 
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VII. UNITARY NOISE AND ERROR CORRECTION CODES 

In this section we are going to study multiunitary noise (6.1), also called random external fields, and look for 
existence of error correction codes, i.e. completely protected subspaccs. In general the number r of unitary operators 
defining the channel can be arbitrary but we will restrict our attention to the cases in which this number is small. 

A. Bi— unitary noise model 

The case in which r = 2, referred to as bi-unitary noise was recently analyzed in [10, 24]. Let us rewrite the 
dynamics in the form 

p' = = qVfV^ + (1 - q)V2pV^ . (7.1) 

and assume that we deal with the system of two qubits. Then both unitary matrices Vi and V2 belong to ?7(4) while 
probability p belongs to [0, 1]. The problem of finding the compression C for the above map is shown to be equivalent 
to the case 

p" = ^*{p)=qp+{l-q)UpU^ (7.2) 

where U = V^V2. 

Thus the error correction matrix a of size two defined by eq. (4.7) reads 



Vz(l - 9)A* 1 - q 

where A is solution of the compression problem for U 

CtC/C = A-l. (7.4) 

Thus to find the error correction space for the bi-unitary model it is sufficient to solve the compression equation for 
a single operator U. A solution exists for any unitary U [10], but for simplicity we will consider here the generic case 
if the spectrum of U is not degenerated. Assume that the phases these unimodular numbers Zi, Z4 are ordered and 
that \ipi) denote the corresponding eigenvectors. 

Let A denote the intersection point between two chords of the unit circle, ziz^ and Z2Z4; compare Fig. 8b. This point 
can be represented as a convex combination of each pair of complex eigenvalues, 

A = (1 - a)zi + az3 = (1 - 6)2:2 + bz^ , (7.5) 

where the non-negative weights read 

a = — — — =: sin^ Oi and b = — — — =: sin^ 02 (7.6) 

Zs — Z\ Z4 — Z2 

and determine real phases 9\ and 62- Note similarity with respect to the construction used in the Hermitian case, in 
which (5.5) represents a convex combination of points on the real axis. In an analogy with the reasoning performed 
for a hermitian T we define according to (5.7) an orthonormal pair of vectors \tpi) and \tp2) and define the associated 
isometry C : tpj. Since (U'V'ijV'i) = A = {Uip2\ip2) and (U'V'i|V'2) = = {Uij2\'>pi) then CUC = XI. Therefore 
A belongs to A2{U) as claimed and the range of C provides the error correction code for the bi-unitray noise (7.2) 
acting on a two-qubit system. 

In the case of doubly degenerated spectrum of U the complex number A is equal to the degenerated eigenvalue, 
so its radius, |A|, is equal to unity. In this case the matrix a given in (4.6) represents a pure state, a = a^, so the 
two-dimensional subspace spanned by both eigenvectors corresponding to the degenerated eigenvalues is decoherence 
free. 

Bi-unitary noise model for higher dimensional systems was analyzed in [24]. It was shown in this work that for a 
generic U of size N there exists a code subspace of dimensionality k = int[(A'' + 2)/3]. This result implies that for a 
system of m qubits and a generic U of size N = 2"^ there exists an error correction code supported on m — 2 qubits. 
Furthermore, if N = d"^ and d>3, there exists a code supported on to — 1 quantum systems of size d. 
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B. Tri— unitary noise model 



Consider now a model of noise described by three unitary operations acting on a bipartite, N = 2 x Nb system, 

p' = = qiVipV^ + q2V2pWl + (1 - 91 - q2)V^pV^ . (7.7) 

Performing a unitary rotation in analogy to (7.2) wc obtain an equivalent form 

/ = $*(p) = qip + q2UipUl + (1 - (?i - q2)U2pUl . (7.8) 

The model is thus characterized by two unitary matrices of size iV, namely Ui = V^V2 and U2 = V-^V^. and two 
weights qi and q2 , which we assume to be positive with their sum smaller than unity. 

To find a simplest error correction code for this model one needs to find a two-dimensional subspace, which forms 
a joint solution of three compression problems 

C'<UiC = \u,l 

C^U2C = Ac/,1 , (7.9) 
C^WC = Xwl 

where W = UlU2- Each of the above three problems may be solved using the notion of the higher rank numerical 
range of a unitary matrix. However, for generic unitary matrices Ui and U2 of size 4 the corresponding compression 
subspaces do differ. Thus for a typical choice of the unitary matrices the tri-unitary noise model will not have an 
error correction code, for which it is required that the subspace C solves all three problems simultaneously. 

There exist several examples of two commuting matrices Ui and U2 of size N = 4., such that they possess the same 
solution C of the compression problem. However, to assure that it coincides with the solution of the same problem 
for W = UlU2, we will analyze an exemplary system of size n = 2 x 3. Consider two unitary matrices of a tensor 
product form. 



U2 = Ua®Ub 



(7.10) 



where 




C/a = I e-^" I and C/b = J ) . (7.11) 

Observe that Ui and U2 do commute, so they share the same set of eigenvectors. Assume that the phases satisfy 
a e (7r/2, tt) and ^ € (O, min{a, 2(7r — a)}) . Then the ordered spectra of both matrices read 

C/i = diag{l,e*«,e'",e*("+«),e-*",e'(f-")}, i/a = diagjl, e'«, e"*", e*(«-"), e'", e*("+f) }, (7.12) 

and differ only by the order of the eigenvalues. Both unitary matrices arc represented in Fig. 9 in which z^, i = 1, . . . , 6 
denote the ordered eigenvalues of U\ while = 1, ... ,6 are eigenvectors of this matrix. The same states form 

also the set of eigenvectors of U2, but they correspond to other eigenvalues. Let z'^ denote the ordered eigenvalues of 
U2- Then \ip^) corresponds to Zg = Z5 while |(/?5) corresponds to Z5 = 03. 
The third of the unitaries also has also a tensor product form, 

W = ulU2 = {U\ ® Ub)HUa Ub) ^Ul® I2. (7.13) 
Hence the spectrum of W , denoted by z", consists of three pairs of doubly degenerated eigenvalues, W = 
diag|l,l,e-2^",e-2'",e2»",e2*"|, see Fig. 10. 

Numerical range of rank two for matrices Ui, U2 and W is shown in the pictured as a gray region. Each point 
A G A2(C/i) offers a subspace C2 which forms a solution of the first of three equations (7.9). However, the other two 
equations restrict further constraints for A. 

To construct an error correction code for the tri-unitary noise model we are going to follow the strategy used above 
for solving the compression problem: we split the Hilbert space into a direct sum of two subspaces of size three, and 
try to construct a single state in each subspace. More formally we define the subspace 



C2=^|V'.)(V'd (7-14) 
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where each state is obtained by a coherent superposition of three eigenstates of Ui, 

iV'i) = Val\(pi) + y/E^lfs) + VaE\(p5) .y-^gN 

Since the unitary operators Ui can be expressed as tensor product of diagonal matrices (e.g. U2 = Ua ® Ub), their 
joint set of eigenvectors consits of product pure states only. On the other hand, the states li/ji) and \ip2) are by 
construction entangled. 



Imz 




b) 



Imz 



Rez 



A,(U,) 




I^4> 



Rez 



FIG. 9: Numerical range (gray space): a) A2(l7i); b) A2((72) 



The weights ai are defined as a weights obtained by representing point A by a convex combination of the triples of 
eigenvalues. Since we wish to get a space C being a joint solution of all three equations (7.9), we are going to require 
that the same weights ai can be used to form the compression value A as a combination of both triples of eigenvalues 
for each spectrum. 




= aizi + 03^3 + 05^5 = a2Z2 + a4Z4 + a^ZQ 

— 012^+0323 + 05^5 = 0222 + 04^4 + 06-26 

= oiz" + 03Z3 + 05Z5 = 02Z2 + 0424 + oeZg 



(7.16) 



where Zi, z[ and z'l denote ordered spectra of C/i, \J2 and W ^ respectively. It is now clear that for a generic choice of 
\J\ and U2 (which implies W = C/^ C/2), such a system has no solutions. However, if both diagonal matrices are of the 
special form (7.12), there exists a solution of the problem. The weights a.i satisfy 



Oi = 02 = 1 



1 



■ cos a 



\ 0-3 



05 



O4 



0-6 



and imply the following compression values 




2 — 2 cos a 
1 

2 — 2 cos a 






2 cos a 



(7.17) 



(7.18) 



Due to the symmetry of the problem the latter number Aw is real. 

Substituting the weights (7.17) into (7.15) we get an exphcit form (7.14) of the subspace C. It is now easy to check 
that this subspace satisfies simultaneously all three equations (7.9) with compression values given by (7.18), hence it 
provides a two dimensional error correction code for this noise model. This solution is correct for any unitaries Ui 
and U2 having any set of eigenvectors \(pi), i — 1, . . . , 6 and spectra given by (7.12) and parameterized by phases a 
and ^. 
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The above construction can be generalized for a tri-unitary noise model acting on larger system of size N = 3 x K 
[31]. An error correction code of size K exists in this case, if matrices Ui and U2 have the tensor product form (7.10), 

where Ua = diag|l, e*", e~*"| as before and Ub — diag|l, e*^^, e*^^, . . . , e*^^ |. The code subspace C ~ X^ili 
is then obtained in an analogous way, by representing the Hilbert space as a direct product of K subspaces of dimension 
three each and constructing each state jV'i) as a coherent superposition of three eigenstates oiU\ corresponding to a 
triple of eigenvalues z;, zij^k and zij^^k for / = 1, . . .isT. Note that the code space constructed here for the bipartite 
system does not have the tensor product structure, since it is spanned by entangled states (7.15). 



Z "=Z " 


Imz 






\ 


K \ 


\ l^'l) 1^2) 

v\z"=Z" 
\j '■I '■2 




i 

j 




/] Rez 
^A,(W) 


z"=z" 

1^5) 1^6) 







FIG. 10; Numerical range A.2{W) is represented by a dark triangle 



VIII. CONCLUSIONS 

This paper concerns finite dimensional instruments or Kraus measurements, acting on a quantum system with 
Hilbert space TL. We have proved a version of Heisenberg's Principle, which connects 'darkness' to 'protection' of a 
subspace £ of 7i. 'Darkness' expresses the lack of visibility of the information contained in C from the measurement 
outcome, and 'protection' the degree to which this information remains present in the quantum system. Complete 
darkness corresponds to complete recoverability of information as in error correction codes. 

We have presented examples of darkness and protection: instruments arising from random external fields, arbitrary 
rank 2 channels, and biased permutation channels. Bi-unitary noise models were analyzed recently in regard to their 
error correction properties in [10, 24]. Here we have also considered tri-unitary noise. For a a certain class of tri- 
unitary noise models acting on a 3 x A' quantum system, we have explicitly constructed an error correction code of 
size K. Although this particular noise model might be considered as not very realistic, we tend to believe that the 
technique proposed can be applied to a broader class of quantum systems. 
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